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We consider indirect detection signals of atomic dark matter, with a massive dark photon which 
mixes kinetically with hypercharge. In significant regions of parameter space, dark matter remains 
at least partially ionized today, and dark atom formation can occur efficiently in dense regions, 
such as the centers of galactic halos. The formation of dark atoms is accompanied by emission of a 
dark photon, which can subsequently decay into Standard Model particles. We discuss the expected 
signal strength and compare it to that of annihilating dark matter. As a case study, we explore 
the possibility that dark atom formation can account for the observed 511 keV line and outline the 
relevant parameter space. 


I. INTRODUCTION 

Atomic dark matter arises in models in which dark 
matter (DM) couples to a dark Abelian gauge force, 
and the relic DM abundance is due to a dark particle- 
antiparticle asymmetry (see e.g. m)- It is a minimal 
implementation of the asymmetric DM scenario, which is 
motivated by the similarity of the observed dark and or¬ 
dinary matter abundances [SHE]. Moreover, atomic DM 
is a simple possibility for self-interacting DM. The lat¬ 
ter has emerged as an alternative to the collisionless cold 
DM paradigm that can resolve the discrepancies between 
collisionless cold DM simulations and observations of the 
galactic structure [TuanHsaQ In fact, because in the 
atomic DM scenario the dark particles couple to a light 
force mediator, the DM self-interactions may manifest 
as long-range, with scattering cross-sections which de¬ 
crease with increasing velocity. This feature implies that 
the DM self-interactions may affect significantly the dy¬ 
namics of smaller haloes, while having negligible impact 
on larger scales p51128j . where the collisionless cold DM 
paradigm has been very successful. However, there is also 
plenty of parameter space in which atomic DM behaves 
as collisionless. 

It has been previously observed that asymmetric DM 
with long-range self-interactions may produce detectable 
indirect detection signals due to the formation of DM 
bound states inside haloes today m- Dark bound-state 
formation is accompanied by emission of a dark force 
mediator, which can subsequently decay into Standard 
Model (SM) particles, via a so-called portal interaction. 
In this work, we explore this possibility within the sce¬ 
nario of atomic DM coupled to a light but massive dark 
photon which kinetically mixes with the hypercharge. 

The atomic character of DM is evident in the case 


^ For an overview of the collisionless cold DM challenges at galactic 
and subgalactic scales, see Ref. m- 


of asymmetric DM coupled to a massless dark photon; 
in Ref. [27], it was rigorously shown that DM remains 
atomic even if the dark U(l) force is mildly broken. In 
the case of a massless dark photon, gauge invariance man¬ 
dates that DM must be multicomponent: The dark asym¬ 
metry generation has to occur via gauge-invariant opera¬ 
tors, therefore the net gauge charge carried by the asym¬ 
metric population of a dark species must be compensated 
by an opposite gauge charge carried by an asymmetric 
population of another dark species. The relic abundances 
of the two dark species may (partially) combine into dark 
atoms in the early universe, as well as inside haloes to¬ 
day. This is, of course, analogous to ordinary protons and 
electrons, whose electric charges compensate each other. 

Reference mi examined in detail the case of the dark 
photon acquiring mass via a dark Higgs mechanism in 
the early universe. If the dark photon mass is small, the 
dark phase transition occurs cosmologically at late times, 
thus leaving unaffected the preceding cosmology, includ¬ 
ing the dark asymmetry generation and the dark recom¬ 
bination. It then follows that - as in the case of a massless 
dark photon - DM has to be multicomponent, with the 
dark ions potentially forming dark atoms in the early uni¬ 
verse. The range of dark photon masses for which DM is 
atomic includes much of the parameter space where the 
DM self-interaction is sufficiently strong to affect the dy¬ 
namics of haloes; thus, the multicomponent and atomic 
character of DM may not be neglected when studying 
this effect m- Moreover, the rich composition of DM 
can result in indirect detection signals, some of which we 
explore in this work. In particular, dark atomic bound 
states may form today from the residual ionized fraction 
of DM. The dark photons emitted in these processes - be¬ 
ing massive - may subsequently decay into SM particles 
via their kinetic mixing with hypercharge. 

Related ideas, involving indirect detection signals from 
excitations and de-excitations of bound states of asym¬ 
metric DM in haloes, have been explored in Refs. [38lI41j . 
Indirect signals from the formation of bound states by 
(weakly-interacting) symmetric DM have been discussed 
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in Refs. [1^44| . 

This paper is organised as follows. In S ec.[TTl we briefly 
review the atomic DM scenario. In Sec. |III[ we discuss 
the signal expected from the formation of dark atoms in 
haloes and we compare it to that of annihilating symmet¬ 
ric DM. In Sec. |IV[ we consider the possibility of the 511 
keV line being produced from the decays of dark photons 
emitted during the formation of dark atoms in the Milky 
Way. We conclude in Sec. |V] 


II. ATOMIC DARK MATTER 


We consider fermionic DM charged under a dark Ud(1) 
gauge symmetry, under which the SM fermions are un¬ 
charged. We also assume that the low-energy theory de¬ 
scribing the dark sector possesses a global U(l) symmetry 
- a dark baryon number - under which DM transforms, 
and that the relic population of DM is due to an excess of 
dark baryons over dark antibaryons. As discussed in the 
introduction (see m for details), if the Ud(1) symmetry 
is unbroken or mildly broken, gauge invariance implies 
that DM consists of two stable particle species, oppo¬ 
sitely charged under Ud(1); we shall call them the dark 
proton, and the dark electron, e^, and assume they 
are fundamental. We denote their masses by irip and m-e, 
and take their Ud(1) charges to be qp = -1-1 and qe = —1 
respectively. We also take nip ^ me- The low-energy 
physics of interest is governed by the Lagrangian 

C=p^ (ilf) - mp)p^ + {il]) - TOe)e„ 

“ + 2^7 1 ( 1 ) 

where with Ad being the dark- 

photon field. (As in QED, we use for the field in the 
Lagrangian, and for the dark photon when discussing 
processes in which it participates.) is the dark pho¬ 
ton mass, which may be generated either via the Higgs 
mechanism or the Stiickelberg mechanism. The covari¬ 
ant derivative for p^ and is + iqjgA^, where 

qj is the respective charge and g is the gauge coupling of 
the dark force. In the following, we more commonly use 
the dark fine structure constant = g^ 

Atomic bound states of dark protons and dark elec¬ 
trons exist if My < where is the p^ — 

reduced mass. 


WpWe 
nip + nie 


( 2 ) 


The ground state of the dark Hydrogen atom has 
mass 


TOjj = nip + ni^ — A , 

where A is the ground-state binding energy, 

1-^V- 



( 3 ) 

( 4 ) 


Note that p^ and satisfy the consistency relation 

4Md < mjj + a , (5) 

where the equality is realized for nip = nie. (In most 
cases, it suffices to approximate ([^ with Ap^ < m^.) 

Provided that My < A, dark atoms can form radia- 
tively, with emission of a dark photon, 

Pn ~b ^ ^n ~b Tr) ■ (^) 

Evidently, the upper limit on My for the radiative forma¬ 
tion of dark atoms is stronger than that for the existence 
of dark atomic bound states. Dark protons and dark elec¬ 
trons bind partially into dark Hydrogen atoms during the 
epoch of dark recombination in the early universe, which 
has been studied in detail in Ref. [53] . The phenomenol¬ 
ogy of DM today is largely determined by the residual 
ionization fraction. 


X 


D 


rip 

Uh + Up ’ 


( 7 ) 


where rij denotes the number density of the j dark el¬ 
ement today. Obviously, rip = rig. In the following, we 
will often denote n^+Up = udm- The residual ionization 
fraction can be approximated by |25j 


min 


1, 10-1° 


Cdr 



/ ni^p„ 1 \ 

^GeV^ PSj 


( 8 ) 


In Eq. ([^ and in the following, the parameter ^ = 
min[l, Tjj/Ty] is determined by the ratio of the dark sec¬ 
tor temperature To, to the ordinary sector temperature 
Ty. ^ may in general vary with time, albeit typically 
mildly. The subscript denotes the relevant time, with 
“DR” referring to the epoch of dark recombination. In 
Eq. <§ we have inserted a phase-space suppression factor 
relevant when A ~ My, 


PS = 


Jl- 

V {/^ + p,vlJ2Y^ 


( 9 ) 


where Urei is the average relative velocity of p^ and 
at the relevant time (for Eq. Q, this is the time of dark 
recombination). Although the approximation of Eq. ([^ 
for the ionization fraction works well when = 1 and 
Xq 1, it is less satisfactory when x^ < 1. However, as 
the ionization fraction depends strongly on , this is a 
relatively small region of parameter space. 

The efficient annihilation of the dark antiparticles {p^ 
and Cj^) in the early universe, sets a lower limit on the 
DM annihilation cross-section, and therefore on a^. In 
the presence of a particle-antiparticle asymmetry, the 
dark antiparticles are diminished to less than 10% (1%) 
of the DM density, if the annihilation cross-section is 
only about 1.4 (2.4) times or more the value required 
for symmetric thermal relic DM [45]. Eor DM consisting 
of Dirac fermions coupled to a dark photon, the cou¬ 
pling required to obtain the observed DM density in the 
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symmetric regime has been calculated in Ref. |46j . tak¬ 
ing into account both the Sommerfeld enhancement of 
the DM annihilation, as well as the formation and de¬ 
cay of dark particle-antiparticle bound states, which also 
contributes to the overall DM destruction rate. More¬ 
over, in the atomic DM scenario, can annihilate 

either into or e^e^, with equal cross-sections; this 

doubles the total annihilation rate. Thus, as the mini¬ 
mum value of a^y for efficient annihilation, it suffices to 
consider that determined for symmetric DM in the pres¬ 
ence of the ^yy^yy annihilutton channel only. We require 

^ a_D,syin(mp), with gS] 

Q^D,sym(^p) / \ 

0.031 ^ 

where the subscript “FO” refers to the time of the DM 
freeze-out. Since me ^ Wp, this obviously also ensures 
efficient annihilation of Cyy. In the next section, we shall 
see that this minimum value of implies both upper 
and lower bounds on the signal strength expected from 
the radiative formation of dark atoms today. 



III. SIGNAL PRODUCTION 
A. Dark atom formation 


The cross-section times relative velocity, for the radia¬ 
tive formation of dark atoms shown in ([^ is m 


((TUreOesF 



c' 

( 1 + 0 ^ 


g—4c;arccot((^) 
1 - 


■PS, (11) 


where C is the ratio of the Bohr momentum to the initial 
state center of mass momentum. 


c=^ = ^, ( 12 ) 

fiVipel Vrel 

and again we have introduced a phase-space suppression 
factor. For the indirect detection signals produced from 
dark atom formation inside haloes, we will be mostly 
interested in the regime Urei OLyy. In this limit, 


, , cm^ 

(crUreOesF ~ 3.5 X 10 

f ayy x3 / lOGeV 

V 1 O- 2 J 



(13) 


occur due to the nonzero mediator mass. Such resonances 
affect only a small region of parameter space, and serve 
to enhance the signal. 

Being Sommerefeld enhanced, the dark atom forma¬ 
tion becomes very efficient in the non-relativistic envi¬ 
ronment of DM haloes. However, the expected signal to¬ 
day depends on the interplay between the strength of the 
interaction cross-section given above, and the early uni¬ 
verse cosmology, which determines the residual ionization 
fraction of DM. We discuss this interplay in Sec. |III C| 

Partial wave unitarity sets an upper limit on the total 
inelastic cross-section; in the non-relativistic regime, this 
is CTinei.j < (2j -f l) 7 r/(/i^u 2 gj) gS] . This gives a rough 
estimate of the range of validity of Eq. (111 ^ 


a. 


< 


OiD,, 


0.5. 


(14) 


B. Dark photon decay 

We have noted that the SM particles are uncharged 
under the Ud(1) symmetry, while the dark fermions are 
uncharged under the SM gauge group. However, the 
dark sector and the SM sector may couple through ki¬ 
netic mixing of the U(l) gauge bosons. We introduce the 
renormalizable operator gilsn] 

/Imix = (15) 

Once SUl(2) is broken, this term induces a coupling be¬ 
tween the dark photon and the SM photon, as well as 
between the dark photon and the Z boson. 

The dark photon may decay into charged fermions, at 
the rate [51] 


= /em • -^chemM^ , 


(16) 


where ctem — 1/137 is the electromagnetic fine-structure 
constant, and /em accounts for the number of kinemati¬ 
cally allowed decay channels. The dark photon may also 
decay into Dv via its mixing with the Z boson, with rate 




3M4 


3 ^ 4cos^(0vv)w^ 


(17) 


Although the decay jyy 77 is forbidden, the dark pho¬ 
ton can decay to three photons via a charged fermion 
loop, at the rate g^j 


which exhibits the familiar oc 1/urei scaling of 
Sommerfeld-enhanced processes at low velocities. This 
low-velocity enhancement saturates when the momen¬ 
tum transfer between the incoming particles becomes 
less than the mediator mass, i.e. at Urei ^ M-y/p-o'i for 
a massive dark photon, in Eq. (13), we may thus replace 

Urel -t max[Urel,M.y/^^]. 

Although we account for the saturation of the cross 
section, we do not consider here resonances which may 


f"7D->'37 ~ 


17e^aE[^ikf® 
2'^3®537r2m® ’ 


(18) 


3 A more precise determination of requires partial wave 

decomposition of the differential cross-section leading to Eq. 
which is beyond the scope of the present work. Such precision is 
not significant for our purposes. 
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where We = 511 keV is the ordinary electron mass. The 
decay to charged fermions dominates unless it is kine¬ 
matically forbidden. 

The cosmology of relic dark photons was considered 
in m- Provided that the UD(l)-breaking phase transi¬ 
tion occurs before BBN, dark photons with a mass above 
1.022 MeV decay before BBN if 


e > 


10-10 

„l/2 

Iem 


10 MeVV^^ 

~W~ ) 


(19) 


Otherwise, we must require that during BBN ^bbn = 
Td/Tv < 0.6, to satisfy constraints on relativistic energy 
density. In this case, further constraints (such as avoiding 
altering the time of matter-radiation equality) may be 
relevant. Experimental bounds on e are explored in IS3 
[54] . For dark photon masses around ^10 MeV, e < 10“® 
due to supernova 1987a, while for dark photon masses 
on the GeV scale, e < 10“^ is constrained by collider 
experiments. 

Therefore, the relic dark photons will typically have de¬ 
cayed before BBN, and for an even wider range of mix¬ 
ings e, before the present day. The indirect detection 
signals relevant today must arise from ongoing produc¬ 
tion of dark photons. Consequently, we focus on dark 
photon emission which accompanies ongoing dark atom 
formation in DM haloes. 


C. Signal strength 


Bound states are produced at the rate per volume 


f \ ^^(^^reOBSF 2 /'on^ 

= (0-llrel)BSF«P»^e = -, (20) 


dV 




where is the DM density, and we took into account 
that 


PDM = ™P + IT-e nie 

= (1 - x^)nr)Mm^ + udm (wp + me) 
~ (1 - x^)nuMm^ + x^ nuM m^ 

= udm m^ , 


with riDM = + iip) s-iid fh® binding energy was ne¬ 
glected. This allowed us to substitute UpUe = = 

^IpUK- 

Evidently, the signals produced by dark atom forma¬ 
tion scale with the DM density in the same way as DM 


annihilation. We define the signal strength s^gp, as[^ 


Sbsf 


(g-UreOBSF 

m? 

H 


( 21 ) 


which facilitates the comparison with signals expected 
from annihilating DM, and contains the combination of 
parameters that is generally constrained by observations. 
In the regime ^ Urei, 


2^TT^a^ X PS 3-10-1® GeV^^^ 

’ a^v^eixPS 

( 22 ) 


where inside the square brackets, the factor on the left 
corresponds to the fully ionized case and the factor on 
right corresponds to the partially ionized regime. The 
energy available to the dark photons emitted in the for¬ 
mation of dark atoms is the sum of the binding energy 
and the kinetic energy of the reduced system, 

(23) 

In most cases of interest, Urei ^ ~ A. 

A contour plot of the signal strength, s^gp, is shown 
in Fig. [2 In the fully ionized regime, the signal strength 
increases with the coupling ; however, in the partially 
ionized regime, the signal strength decreases rapidly with 
increasing . This is a consequence of the dramatic 
decrease in the ionization fraction. The sharp line at 
cCp, = 1 is a feature of the approximation of Eq. Q ; an 
exact ionization fraction would round the sharp corners. 
The consistency relation of Eq. implies that at every 
point in the {a ^, m^fi ^) plane, there is a maximum value 
of the energy imparted into the dark photon which is 
attained for nip = nie and can be expressed as 

Wy < ^ («„ + I'rel) ■ (24) 

Since the radiative capture of to bound states is 

possible only for this condition constrains also 

the possible dark photon decay channels and the indirect 
detection signatures. In Fig. the orange dashed line 
marks = 1.022 MeV; to the left of this line, the de¬ 
cay 7 ^ —)■ e"^e“ is kinematically forbidden for any choice 
of My < ccy^o- (However, the dark photon may still de¬ 
cay to neutrinos or three SM photons.) We see that if 


^ More properly, 5ggp should be defined in terms of the average 
(crUj.ei)BSF over the velocity distributions of the participating 
particles. For Urei "C , which is the regime of greatest in¬ 
terest, ((ji^reOBSP oc l/'T’reh for a Maxwellian distribution, 
(l/urei) = y/6/7r/firei ~ l.SS/fireij where firei is the rms value 
of the relative velocity. However, for simplicity, we will use 
(l/urei) 1/^reh denote firei with Ufei• A proper averaging 
would enhance the expected signal. 
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upper limit on given in (25), bound the expected sig¬ 
nal strength from both above and below. 


For convenience, we first invert Eq. (101 numerically, 
to obtain the maximum dark proton mass for which anni¬ 
hilation is efficient in the early universe, given a coupling 
®D) ^ Wp_sym(oD)- We find the approximation 


Tod 


) 


TeV 


o.osUf^" 


1.7 


-I 0.42 


1 + 


Ctr 


o.osUp/' 


(27) 

Furthermore, since TOp_sym(ajD) ^ Wp ^ = 4A/q;^, 

we may employ the above to obtain a lower bound on 
as a function of the binding energy, > aD_sym(A). 
We find numerically 


bD,sym(A) 


240 MeV 


1 -k 


( 

I 240 MeV 


0.5 


n -0.13 


(28) 


FIG. 1. Countour plot of the signal strength s^gp, from bound 
state formation, in units of GeV”"^. (For convenience, we note 
that GeV”"* ~ 10“^^cm^s“^ TeV“^.) We have used = 1, 
Urei = 200 km/s. We assume that the dark photon mass is 
less than the binding energy, < A, such that the radiative 
formation of bound states is possible. To the left of the orange 
dashed line, the decay —>■ e^e“ is forbidden, because the 
binding energy is A « a^^p,/2 < 1.022 MeV for all possible 
values of Above the orange dashed line, whether the 
decay —>■ e''"e“ is kinematically allowed depends on M.y. In 
the grey shaded region, the coupling is insufficient to remove 
the symmetric DM component in the early universe, for any 
consistent choice of (Eq. (25l). 


the dark photon decay into charged fermions is possi¬ 
ble, the greatest signal strength that can be produced 
is Spgp 10“® GeV“^. More generally, while the signal 
strength depends on various parameters and is not fully 
determined at a specific value of the dissipated energy, 
it is bounded both from above and below. We shall now 
explore the expected signal range. 


D. Circumscribing the signal strength 


Equations ( ^ and ( |14[ ) suggest that the DM coupling 
to the dark pt 


rotons should be in the range 

C^i^jSym(^Up) < Op, < • 


Moreover, 


(25) 


8A/al = 4^ 


< 


m. 


5, 2 mn 


(26) 


We shall now employ Eqs. (25) and (26) to circumscribe 


the strength of the signal expected from dark atom for¬ 
mation in haloes today, given by Eq. (21); in Sec. HIE 


we compare it to the signal expected from annihilating 
symmetric DM. As we shall see, both the lower and the 


1. Signal range at fixed 


From Eq. (22), it follows that for a fixed value of Op 


there is a maximum signal strength that can be produced, 
which is given by 


SBSF,max(ap,) 


10-19 QeV-4 ^2^ 
a® Urel 


(29) 


provided that » 'Crei and < A. (For lower val¬ 
ues of Op,, an expression for seSF.max can be obtained 
using Eq. (111.) This maximum signal is produced if 
Coj,(mHAio)7GeV^ < lO^o and is shown as a function 
of Op, by the blue solid line in Fig. 

For larger values of m^/ip,, the signal strength de¬ 
creases; however, from Eqs. ( [^ and (27) we find, 
'a < ml < ml^^^a^). This implies a 


mr-u ^ rnHA < ml < 
lower bound on s^gp. 


'^BSF,niin (oip,) 


2977^ al 

77i 4 ^y^(a^) Vrel ' 


(30) 


In Fig. we sketch the signal strength s^gp as a func¬ 
tion of Op, for different values of m-^ji^. In the orange 
region, the decay y^ —>■ e+e” is forbidden for all choices 
of /ip, and TTip, which satisfy the consistency relation ([^. 
In Fig. we plot both the maximum and the minimum 
signal strengths expected from bound state formation, 
over a larger parameter region. We emphasize that for 
fixed Op, the signal strength depends only on ^dr; how¬ 
ever, the minimum signal depends on ^ts minimum 
coupling Op, which permits sufficient annihilation of the 
symmetric component does. 


2. Signal range at fixed DM mass, rUp 


The upper bound on s^gp corresponding to the par¬ 
tially ionized branch of Eq. (21), can be re-expressed in 
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FIG. 2. Maximum signal strength s^sp, for = 1, 

Wrei = 200 km/s. We assume that the dark photon mass 
is always lower than the energy dissipated in the forma¬ 
tion of bound states, such that the latter can occur radia- 
tively. In the orange region, the decay ^ e''"e“ is for¬ 

bidden for all of parameter space. The dashed lines corre¬ 
spond to the signal in the fully ionized region. (Note that 
GeV-'‘ ~ 10-“cm3 s"^ TeV^) 


terms of m^, using Eq. (10) and irip ^ 


3 • 10-GeV-4 eL 


•^BSF.max ) 


(31) 


Moreover, using the unitarity limit of Eq. (14), and also 
the consistency condition © , we obtain a lower limit on 


Sggp as a function of , 


■^BSF.min (tU-jj ) ~ 

. 3 • 10-19 GeV-^CL 

iin - — - - 

3e-^(m^/4) Vrel ’ afj.uni^^rel 


(32) 


We sketch Eqs. (31) and (32) in Fig.[^ 


3. Signal range at fixed value of the dissipated energy 

The detectability of signals from DM related processes 
depends of course on the energy of the relativistic prod¬ 
ucts of these processes. In contrast to the more familiar 
case of DM annihilation, the energy dissipated during 
dark atom formation is only a fraction of the dark parti¬ 
cle mass, ~ A ^ . We shall thus now express the 

minimum and maximum signal strengths expected from 
dark atom formation in terms of the binding energy, A. 


FIG. 3. Minimum (red dotted) and maximum (blue dashed) 
signal strength from dark atom formation, as a function of the 
dark hne structure constant, . The green solid line is the 
expected signal strength from the annihilation of symmetric 
DM whose coupling strength and mass are related by the 
observed DM density (see text for details). We have used 
= ^pQ = 1 (which holds if the dark sector temperature is 
greater than the visible sector temperature at both times) and 
Urei = 200 km/s (smaller Urei would enhance s^gp). For atomic 
DM, we assume that the dark photon mass is 0 < M.y < A, 
such that the radiative formation of bound states is possible 
and the decay of dark photons is kinematically allowed. (Note 
that GeV-'‘ ~ 10-“cm®s-^ TeV-^) 



FIG. 4. Same as in Fig.[^ for the signal strength as a function 
of the DM mass. For atomic DM, is the dark Hydrogen 
mass. 
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It will be convenient to rewrite Eq. (22), exchanging 
for A, as follows 


2Va7 3-10-19 GeV'^e^ 


3e^ ’ 


a^^rel 


For a fixed value of m^A, the maximal signal is 

r 2 N 5/6 


'5BSF,max 

and occurs if 


3.6 X 10-11 GeV-^i ^ ^ 


/ m^A \ 

VGeVV 


Urel 


5 X 109 a® . 


(33) 

(34) 

(35) 


From Eq. ( [^ , > SA/m^, which together with 

Eq. (35) implies CD4^(’7^H/G^eV) > 40(A/MeV)i/9 along 


the line of maximal signal. Equation (|34|) now becomes 



SBSF.max « 5.3 X 10-1° GeV 


^rel \ ^ / 


This is valid provided that SA/m^ < ajj ^JJ;. For 


> 8A/c 


D.uni’ 


the fully ionized branch of Eq. (33) can 


give a stronger bound. Gollecting everything, we find 


SBSF.max(A) « min 


275^ g 

3e4 Urel A4 ’ 


5.3 X 10 


-10 GeV-^ 


11/8 

DR 


/MeV\ 


5/4' 


l^rel V A / 


(37) 


We obtain lower bounds on s^gp by considering the fully 


and partially ionized branches of Eq. (|33|) separately, and 
taking into account that m 
obtain 


< 2y7ip ^ 2?Up syip, (Hp, ) . AVe 


FIG. 5. Same as in Fig. for the signal strength as a 
function of the energy dissipated into radiation. For bound 
state formation, the dissipated energy is approximately equal 
to the binding energy 5E ~ A, while for DM annihila¬ 
tion 5E ~ Note that thermal relic DM with mass 

tUdm /S MeV, necessitates ^ < 1, otherwise it would 
become non-relativistic and freeze-out around or after BBN 
(dot-dashed green line). 


symmetric thermal relic DM0the latter includes the case 
of WIMP DM, however symmetric thermal relic DM may 
also reside in a hidden sector [^. The rate of annihi¬ 
lations per volume is dTann/dE = Sann, where for 
non-self-conjugate DM 


_ {(TVrel) ann 

4m9 


(39) 


SBSF.min(A) « min ■ 


3 • 10-GeV-^F 


-4^2 


0,1 


1 '^rel 


3e^ Vre\ 


m. 


p,sym 


(«d). 


(38) 


.(A). 


We sketch Eqs. (37) and (38) in Fig.[^ 

We note that Eqs. (29) - (32), ( [^ and (38) assume 
ttp, > Urei; for < Vrei, similar considerations can be 
carried out, albeit the signals from bound-state formation 
weaken significantly. In the above, we used the same 
symbols (sBSF.min and sesF.max) for different functions 
to ease the notation. 


E. Comparison to annihilating dark matter 

We now compare the expected signal from dark atom 
formation, with the expected signal from annihilation of 


For s-wave annihilation, the cross-section which yields 
the observed DM density is [59] 

(o-Urel)ann ^p^ X 4.4 X IQ-^® Cm®/s . (40) 

Here we consider only DM annihilating via a short-range 
interaction; for DM interacting via a long-range force, the 
annihilation cross-section is enhanced at low velocities by 
the Sommerfeld effect, which may affect both the freeze- 
out [milB] and the annihilation signals today [4211441160] . 
We compare Sann to SBSF.min and SBSF.max in Figs.[^~[^ 


Asymmetric DM may also result in suppressed but detectable 
annihilation signals due to the subdominant population of an¬ 
tiparticles left over from the early universe |55l I56| . Moreover, 
annihilation signals are predicted in scenarios where DM pos¬ 
sessed a particle-antiparticle asymmetry at early times that was 
subsequently erased due to oscillations EZ]. We do not consider 
these possibilities here. 

































In Fig. we present SBSF.min, SSSF.max and Sann vs. 
the DM coupling which is responsible for the process that 
produces the expected signals. For atomic DM, this is the 
dark fine structure constant . To facilitate the com¬ 
parison, for the symmetric DM case, we shall consider 
Dirac fermions coupled to a massive dark photon that 
is lighter than the DM particles albeit sufficiently heavy 
such that the Sommerfeld effect is not relevant. In this 
case, (crurei)ann ~ Then, the required value 

of the annihilation cross-section quoted in Eq. (401 im¬ 
plies a relation between and which we use in 

sketching the Sann curve in Fig. We see that, for the 
same coupling the signal expected from the forma¬ 
tion of dark atoms exceeds significantly the signal ex¬ 
pected from the annihilation of symmetric thermal-relic 
DM. This is partly due to the Sommerfeld enhancement 
of bound state formation, which renders it very efficient 
at the low relative velocities inside DM haloes. 


Similarly, when compared at equal values of the DM 
mass, dark atom formation appears to be more efficient 
than the annihilation of symmetric DM in much of the 
parameter space, as seen in Fig. This is in part be¬ 
cause for symmetric thermal-relic DM, the DM coupling 
strength is fixed by Eq. (40); in contrast, in the asymmet¬ 
ric DM scenario, stronger annihilation does not diminish 
the DM abundance and there is no upper bound (but 
rather only a lower one) on the DM coupling from this 
consideration. Conversely, a given value of may cor¬ 
respond to much smaller DM mass in the case of asym¬ 
metric DM, than in the case of symmetric DM; smaller 
mass implie^in turn, larger interaction cross-section, as 
seen in Fig. ^ Moreover, cr^gp oc while iTann oc 
and in general, may be much smaller than m^. 


In Fig.[^ we compare the signal strengths at equal val¬ 
ues of the dissipated energy. For atomic DM, 5E ~ A <C 
mjj, while for annihilating DM, 5E = At small 

e, the annihilation of symmetric thermal relic DM is ex¬ 
pected to give stronger signals than dark atom formation, 
despite the latter being a more efficient process. This is 
because dark atom formation is suppressed by the small 
number density of the dark particles, due to the much 
larger DM mass (for the same value of 5E). However, 
annihilating species cannot be cold thermal relics if their 
mass is significantly below ~ 10 MeV, unless ^ <C 1. 
It is thus unlikely that the annihilation of thermal relic 
DM can explain sharp features in the low-energy part of 
the 7 -ray spectrum, such as the observed 511 keV and 
3.55 keV lines. On the other hand, thermal relic DM 
can potentially account for such features if it is asym¬ 
metric and possesses internal structure which gives rise 
to level transitions pHIl40| . such as formation of bound 
states. We discuss in particular, the possibility of the 
511 keV line being produced due to the radiative forma¬ 
tion of dark atoms, in the next section. At larger values of 
the dissipated energy 6E (which also implies larger DM 
mass and larger coupling), the signals from dark atom 
formation can be potentially comparable (or stronger) to 
those expected from symmetric DM annihilation as seen 


in Fig. 

The above discussion does not capture, of course, the 
entire complexity of the indirect detection of DM. The 
detectability of the radiation produced by DM-related 
processes depends not only on the rate and the total en¬ 
ergy dissipated, but also on the energy spectrum and 
the nature of the end products. In the radiative forma¬ 
tion of dark atoms shown in ([^, after averaging over 
the polarizations of the massive photon, the dark photon 
decay products are equally likely to be emitted in any 
direction with a flat energy spectrum, extending to en¬ 
ergies A ~ A. In contrast, the direct products of DM 
annihilation have a sharply peaked spectrum at energy 
E ^ However, most direct WIMP annihilation 

products, e.g. bb, are followed by cascade decays which 
result in an extended final spectrum at lower energies. 
In addition, propagation effects can substantially mod¬ 
ify the spectrum shape for both WIMP annihilation and 
dark atom formation. 


IV. THE 511 KEV LINE FROM DARK ATOM 
FORMATION IN THE GALACTIC CENTER 


In this section we focus on charged decay products of 
the dark photon, which are more easily detectable than 
neutrinos. (The 7 ^^ —>■ 87 decay is extremely suppressed, 
as seen from Eq. (18).) We saw above that, when the 
dark photon decay into charged particles is kinematically 
allowed, at A > MeV, the largest signal strength occurs 
near ^ 0.015 and ^ 400 GeV^. The con¬ 

sistency condition then implies < 10 GeV and 
correspondingly, 40 GeV. (Of course, allowing for 

somewhat lower signals broadens the parameter range.) 
The requirement 1.022 MeV < < A ensures that 

close to the maximal signal region, appreciable signals are 
produced only in the e“'"e“ channel. We shall now explore 
the possibility that in this region - of near-maximal sig¬ 
nal and low positron injection energy Einj < A/2 ~ MeV 
- the formation of dark atoms may explain the 511 keV 
line observed in the center of the Milky Way m- 


A. The 511 keV line 

The 511 keV photon flux from the bulge of 
the Milky Way is observed to be (1.05 ± 0.06) x 
10 “^ photons cm“^ s“^, with a spatial extension of ~ 8 ° 
(FWHM) [61]. This corresponds to an annihilation rate 
of non-relativistic positrons of 

Te+e-.obs^ 1-5 X 10^3 s-S (41) 

within radius r^ax ~ (8°/360°)7ri?sc — 0.6 kpc from the 
galactic center, where i?sc — 8.5 kpc is the radius of the 
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solar circle. E] 

Positrons injected in the Galaxy at some higher en¬ 
ergy, propagate in the interstellar medium, lose energy, 
and annihilate with electrons. The majority of positrons 
survive until they become non-relativistic and annihilate 
at rest to produce the 511 keV line. However, a portion 
of the positrons annihilate while still relativistic. The 
survival probability is Psurv — 0.95 [5H[B3], with a very 
mild energy dependence which we shall ignore. 

The in-flight positron annihilations can produce a sig¬ 
nificant flux of 7 -rays; the observed y-ray continuum 
then constrains the injection energy of the positrons re¬ 
sponsible for the observed 511 keV line. Monoener- 
getic positrons should be injected at energies < 
3 MeV [52], although this bound could be relaxed if un¬ 
certainty greater than 30% in the diffuse y-ray flux is 
incorporatedE] In addition, a broader positron injection 
distribution, such as that produced by the dark photon 
decay, could extend to somewhat higher energies. 


B. Positrons from dark atom formation 


While astrophysical explanations for the production of 
the galactic positrons have been proposed [55H55] . the 
origin of the 511 keV line remains a mystery. Dark mat¬ 
ter related explanations have also been put forward (see 
e.g. [55]). Here we show that the decay of dark photons 
emitted in the formation of dark atoms in the center of 
the Milky Way, can account for the observed positron 
flux. 

The rate at which dark atom formation can contribute 
to the non-relativistic positron annihilation in the center 
of the galaxy is 

,BSF "^BSF ^ I ^ Esuyv , (42) 

where 

I = JdVpl^, (43) 


with the integration extending to radius from the 
center of the galaxy. We parametrize the DM density 
profile by 


Pdm (^) 



(44) 


where = 20 kpc and po is determined such that 
p{Rsc) = 0.4 GeV/cm^. The Navarro-Frenk-White pro¬ 
file corresponds to a = 1, /3 = 3, and y = 1. However, 


^ In deducing the positron annihilation rate of Eq. from the 
observed photon flux, the positronium fraction in the galaxy was 
taken into account. For details see [61[I62| . 

® A weaker constraint, Ejnj $ 20 MeV, arises from the internal 
bremsstrahlung photons associated with any process which pro¬ 
duces charged particles ED- 


the slope may be steeper in the galactic center due to the 
supermassive black hole [70] and/or baryonic matter [71]- 
iza, as has been suggested by numerical simulations m- 
We fix Of = 1 and P = 3, but consider 1 < y < 1.4. 
For this range of y, / ~ (0.1 — 3.3) x 10®® GeV^ cm“®. 
Then, the signal strength required to produce the rate of 
Eq. ([dT]) is 


1.4x10"® GeV 


,_4 f 10®® GeV^ cm"® 


• (45) 


In Fig. [^ we show the boundaries of the parameter 
space that can account for the observed excess, for var¬ 
ious choices of dark photon mass Mj, and the ratio of 
the dark to visible temperatures For each set of val¬ 
ues (/i^, TOjj), we look for the value of that yields 
the observed rate. Because the signal strength produced 
by dark atom formation is bounded, as discussed in 
Sec. HID there is not always a value of that can pro¬ 
duce the observed rate. In Fig. [^ we also require that 
the consistency relation ([^ is satisfied, and that dark 
atoms can form radiatively, i.e. < A. Note that in 
this region of parameter space, the phase-space suppres¬ 
sion factor is not negligible, and we do incorporate it in 
our calculations. These plots do not include astrophysi¬ 
cal and cosmological constraints, which will be discussed 
below. 


Since dark atom formation is more rapid in more con¬ 
tracted profiles (large y), the parameter space which can 
account for the 511 keV flux is correspondingly larger. 
Moreover, a larger temperature ratio Tu/Ty, increases 
the residual DM ionization fraction and strengthens the 
expected signal, thus proving larger parameter space for 
explaining the observed flux. The available parameter 
space is also dependent on the mass of the dark photon; 
if Mj is close to the binding energy, the phase space sup¬ 
pression of the signal limits the parameter space which 
can account for the line. 


For a fixed reduced mass p,^ and bound state mass m -^, 
there are typically two couplings which can produce 
the observed 511 keV flux; one corresponds to fully ion¬ 
ized DM today, while the other corresponds to partially 
ionized DM today. The couplings necessary to produce 
the observed flux are generally smaller in the partially 
ionized branch; therefore, less parameter space is elimi¬ 
nated by the constriant A > and consequently, the 
parameter space in which the observed 511 keV excess 
can be produced is generally larger for the partially ion¬ 
ized branch. 


We must also ensure that the dark photons emitted 
during the formation of dark atoms in the galactic cen¬ 
ter, decay into e^'e" within r < rmax- From Eq. (nq, 
the dark photon decay length is A ~ / {3/^)^ , 

where we took into account the dark photon Lorentz 
boost A/M^. Requiring conservatively A < Fmax/lO — 
0.06 kpc, sets a very comfortable lower bound on the 
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C. Astrophysical and Cosmological Constraints 

Having outlined the region of parameter space in which 
the observed 511 keV line can be produced, we now con¬ 
sider the relevant astrophysical and cosmological con¬ 
straints, and present them in Figs. ^and[^ for the fully 
ionized branch, and Figs. [ 8 | and [TO for the partially ion¬ 
ized branch. 

The most stringent constraint comes from CMB data. 
Bound state formation during recombination deposits ad¬ 
ditional energy into the CMB, which may cause distor¬ 
tions m- The rate of energy deposition from bound 
state formation is 

+ (47) 

where pc is the critical density of the universe, ~ 

0.26, z is the redshift, and 

Pbsf( 2 ) = /(^)Sbsf (48) 



FIG. 6 . Parameter space which can produce the 511 keV 
line, for an appropriate choice of which depends on m^, 
p,^. The dark photon mass is — 2 MeV (top) and M-^ = 
9 MeV (bottom). The solid lines correspond to fully ionized 
DM, while the dashed lines correspond to partially ionized 
DM. Note that the dashed lines encompass all of the region 
enclosed by the solid lines, for the same ^ and 7 , albeit the 
values of are different in the two cases. (As seen from 
Fig. 0 the maximum value of which can produce a 

given signal Sbsf> 1 ® same for fully and partially ionized 
DM.) More contracted halo profiles (larger 7 ) and larger dark- 
to-ordinary temperature ratios ^ in the early universe, lead to 
a stronger signal and hence a larger parameter space. 


kinetic mixing 




( 46 ) 


Here, f{z) is a (redshift-dependent) factor which de¬ 
scribes the fraction of energy absorbed by the CMB; for 
the two-stage process of dark photon emission and de¬ 
cay into e+e“, f{z) ~ 0.45 [T^- Note that in Eq. (48) 


we have used , which incorporates the (possibly non- 
maximal) residual ionization fraction of DM, since dark 
recombination is expected to precede ordinary recombi¬ 
nation. We evaluate s^gp by replacing Urei -^ 7/^01 
noting that the non-zero dark photon mass imposes a cut¬ 
off on the Sommerfeld enhancement of (crUreOesF! dis¬ 
cussed below Eq. (13). At the time of CMB, the 
relative velocity would be ~ a/3^ / p^, were the 

dark ions in equilibrium with the dark photons; how¬ 
ever, Urei is in fact significantly lower, since the dark ions 
have already decoupled from the dark photons [25]. For 
> 2 MeV, the Sommerfeld enhancement of (o'Ui.eOesF 
has indeed saturated. We require m 


(49) 


The spectrum of the injected positrons from the decay 
of the dark photons emitted in the formation of dark 
atoms, is flat and extends between the energies 


UJry 

T 


1 T ^(1 - M2/cc2)(l - 4m2/M2)J ^ 


(50) 


where ~ A is the energy of the dark photon. We 
require A < 20 MeV, such that the average injection 
energy of the positrons is < 10 MeV. This is consistent 
with the constraints derived from the diffuse 7 -ray flux 
discussed in Sec. |IV A[ particularly considering the non- 
monoenergetic positron injection spectrum. However, in 
the parameter space which can account for the 511 keV 
flux and satisfies the CMB constraints, we find that A is 
even smaller. 

The DM self-scattering inside haloes may affect the 
dynamics of galaxies. In smaller haloes, the effect of the 
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Mj = 5MeV, f = I, 7 = 1.4 


My = 5MeV, f = I, 7 = 1.4 




liD [GeV] 

My = 2MeV, f = 1, 7 = 1.4 


liD [GeV] 

My = 2MeV, f = 1, 7 = 1.4 






FIG. 7. The parameter space which can account for the 
511 keV line (white region), within the regime where DM is 
fully ionized, = 1, and for a contracted NFW profile with 
7 = 1.4; the dark photon mass M^, and the early-universe 
dark-to-ordinary temperature ratio are denoted on the top 
of each plot. We have applied various constraints, as described 
in the text and by the corresponding labels. In addition, the 
blue dashed line encloses the region where DM self-scattering 
can affect the dynamics of dwarf-galaxy-sized haloes. Con¬ 
tours of constant binding energy A are also shown. 


FIG. 8. Same as in Fig.[^ for the regime where DM is partially 
ionized, < 1. 
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M^ = 2MeV, f = 1 , 7 = 1.3 



FIG. 9. Same as in Fig. (fully ionized dark matter), for a 
contracted NFW profile with 7 = 1.3. 


M,, = 2MeV, f = 1 . 7 = 1.3 



FIG. 10. Same as in Fig. (partially ionized dark matter), 
for a contracted NFW profile with 7 = 1.3. 


DM self-interactions can bring theory in better agree¬ 
ment with observations. This parameter space is de¬ 
noted by the blue dashed line, which encloses the re¬ 
gion in which the cross section is greater than 0.5 cm^/ 
g at 10 km/s. Currently, the most stringent constraint 
on the DM self-interactions is thought to arise from the 
observed ellipticity of galaxies of the size of the Milky 
Way [IHl dn] ■ We have shaded brown the region where 
CTscat/wDM > 1 cm^/g at 220 km/s. This follows the 
analysis of m- 

Moreover, we require that the dark force provides effi¬ 
cient annihilation of DM in the early universe, with the 
lower bound on given in Eq. (10). This does depend 
weakly on ^po- We note that the decoupling of the dark 
and visible sectors depends on the kinetic mixing e and 
dark electron mass; however, if the two sectors decoupled 
kinetically before DM freeze-out, then ^po will not gen¬ 


erally equal ^dr- However, ^ is slowly varying and the ^ 
dependence of this constraint is rather mild. Therefore, 
for simplicity, we have taken ^po = Cdr- 

In the entire parameter space of interest, is well 
below the unitarity bound of Eq. ( [T^ . We enforce the 
consistency condition of Eq. ([^, and we exclude the re¬ 
gions where > A, since dark atoms could not then 
form via emission of a dark photon. 

Successful BBN constrains the radiation present in the 
universe at temperatures Ty = Tbbn ~ 1 MeV. A rela¬ 
tivistic dark plasma consisting of dark photons, present 
at BBN, would need to be at temperature Td < 0.6Ty. 
Dark photons with mass > 2 MeV would be quasi- 
relativistic or non-relativistic at BBN and this constraint 
(C ^ 0-6) is at least somewhat relaxed. Nevertheless, to 
retain ^ < 1, it is necessary to ensure that the dark and 
the ordinary sectors do not equilibrate via the kinetic 
mixing e. This poses the condition < 10’ -19 

where we assumed that the lightest charged particle in 
the dark sector is the dark Higgs which gives mass to 
the dark photon, and has itself only a somewhat larger 
mass than the dark photon. (If the dark photon acquires 
its mass via the Stiickelberg mechanism, then is the 
lightest dark charged particle, and the above constraint is 
significantly relaxed.) Direct detection experiments sug¬ 
gest a similar upper bound on e. (See e.g. Ref. [7S1 ITT] , 
although the bounds are somewhat strengthened by more 
recent experimental data m- However, a dedicated 
analysis of direct detection data for atomic DM, which 
would take into account the different DM elements and 
the various possible DM-nucleon interactions, is needed.) 
These upper bounds on e are very comfortably compat¬ 
ible with the lower bound estimated in Eq. ( |46[ ), which 
ensures prompt decay of the dark photons inside the halo. 


V. CONCLUSION 

The asymmetric DM scenario allows for rich dark sec¬ 
tor microphysics, which can yield distinct detection sig¬ 
natures. Because in this scenario, the DM relic abun¬ 
dance is protected by a conserved particle-number excess, 
the observed DM density does not set an upper bound on 
the DM couplings to lighter species, despite DM being 
a thermal relic. Dark matter may then possess signifi¬ 
cant couplings to light force mediators; this is, in fact, 
motivated by the observed galactic structure, which cur¬ 
rently can be explained better by DM with sizable self¬ 
interactions rather than by collisionless cold DM. Sim¬ 
ilarly to ordinary matter, DM with long-range interac¬ 
tions may possess rich phenomenology. The formation 
of bound states in particular, is an important feature of 
many such theories. 

In this paper, we have explored indirect detection sig¬ 
nals arising from the formation of atomic bound states in 
haloes today, by asymmetric DM which couples to a light 
dark photon. Level transitions, such as bound-state for- 
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mation, can yield low-energy signals which cannot oth¬ 
erwise be easily produced by annihilating thermal-relic 
DM. However, level transitions may also produce signals 
in the spectrum expected from WIMP annihilation. In 
our analysis, we explored the expected signal strength 
from dark atom formation generically. We also showed 
that the radiative formation of dark atoms can account 
for the observed 511 keV line, provided that the DM pro¬ 
file in the central kpc of our galaxy is steep. The observed 
radiation backgrounds can constrain the parameter space 
of atomic DM with kinetic mixing to hypercharge, based 
on the radiation emitted from the formation of bound 
states, and so also allow correlation of direct and indi¬ 
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